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1 Introduction 



Our understanding of the low-energy behaviour of Quantum Chromodynamics calls for 
the developement of non-perturbative methods for quantum field theories. Some succes 
has been obtained by adapting to relativistic quantum fields well-known non-perturbative 
methods used in the nuclear many-body problem. For instance, variational methods using 
Gaussian wave functionals (which are analogous to the Hartree- Fock-Bogoliubov (HFB) 
kind of approximations) have been applied to self-interacting bosonic field theory [1] , and 
to gauge field theory [2]. One important question, which is currently intensely debated in 
the literature, is how do such non-perturbative approximations respect the symmetries of 
the theory. 

In this regard, we would like in the the present work to revisit this important 
question within the framework of yet an other well known and equally successful method 
from the many-body theory, the random phase approximation (RPA) and its extension, 
the renormalized RPA [3, 4]. The first of these two has experienced in the past quite a 
number of applications in various fields either in condensed matter or the nuclear problem. 
The second one on the other hand has, since its first formulation by Rowe [5], attracted 
very little interest till very recently. Beside its use in the conventional nuclear problem, 
there has been indeed several attempts as of lately to extend its application to quantum 
field theoretic models as well [6, 7]. It is worth reminding that in the standard RPA 
approximation, two-body expectation values are evaluated in the HFB mean-field vacuum. 
In the renormalized RPA approximation, however, these two-body expectation values 
are determined using the true RPA vacuum, inferring to the renormalized RPA a self- 
consistent character. 

Recently the quasi-particle RPA approach has been applied to the linear sigma 
model in order to obtain a correct description of the global chiral symmetry broken phase. 
The standard mean-field HFB approximation gives a finite mass to the pion. However, it 
has been shown that in the RPA approximation the pion is massless in accordance with 
Goldstone theorem [8, 9]. A careful study of the finite temperature chiral phase transition 
in this model reveal, however, that the HFB-RPA approximation leads to a first order 
phase transition. This clearly is an artefact of the approximation as it is well admitted 
that the chiral transition in this model is of second order. This problem can be traced back 
to the fact that the finite temperature induced transition does not happen in the RPA 
vacuum, but rather in the self-consistently built ground-state, namely the HFB state. The 
latter is clearly a wrong vacuum for the theory as it doesn't possess a valley in the broken 
phase. The RPA fluctuations which were crucial in correcting for these shortcomings 
[8, 9] are of no use in the present situation since these are implemented perturbatively. 
Therefore, it is obvious that such an approach is bound to fail in describing the theory in 
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the vicinity of the phase transition. A possible solution to this problem might very well be 
in relaxing the quasi-boson approximation inherent in the standard RPA approach. This 
leads to a self-consistent version of the RPA very difficult to put in practice. Therefore, 
it is more convenient to consider rather the renormalized RPA variant. This is what we 
propose to study here. 

The renormalized RPA approximation has been considered in [6] in the context of 
scalar A0 4 field theory in 1 + 1 dimensions. In this work, the formulation of the RPA 
equations was based on the Dyson equation approach. A preliminary work applies the 
renormalized RPA to O(N) field theory [7]. Here, we will focus on the two-dimensional 
quantum anharmonic oscillator with an 0(2) symmetry. This purely quantum mechan- 
ical model allows to obtain analytical expressions without the problems of divergences 
occurring in quantum field theories. It shows the possibility to have a vacuum state with 
spontaneous broken symmetry. It has also the advantage that we can compare the results 
obtained with the standard and renormalized RPA approximations with exact numerical 
results. This numerical investigation shall be published in a future work. Here we will 
concentrate on the formal aspects of the theory. A crucial point is to check the existence 
of a zero excitation frequency above the vacuum state with broken symmetry in our non- 
perturbative approximations ( this corresponds in quantum field theory to the Goldstone 
mode). Since these quantum mechanical systems are used for demonstrational purpose, 
we will disregard all problems related to the infrared divergences occurring due to the 
presence of this zero energy mode. 

The paper is organized as follows. In the first part, we derive the renormalized 
RPA equations using the equation of motion method [5]. In the second part, we use a 
time-dependent approach. In this formalism, the standard RPA frequencies are obtained 
as small oscillations around the static solution in the time-dependent Hartree-Bogoliubov 
equation. Within this second approach, we are not able to go beyond the standard RPA 
approximation. However, it is interesting to draw the parallel between the two approaches. 



2 Renormalized RPA equations from the equation of 
motion method 

The Hamiltonian for the 0(2) anharmonic oscillator reads : 

H = ^ + ^ + \\X\ + (X 2 + < X 2 >f] +g[Xf + (X 2 + < X 2 >ff + V (X 2 + <X 2 >), 

(2.1) 

where we have considered an explicit (rj ^ 0) symmetry breaking and a spontaneous 
(< X 2 >^ 0) symmetry breaking along the X 2 direction. In the \i < case, the potential 



3 



has a "Mexican hat" shape. 

The case of two particles in a harmonic potential and coupled by a linear interaction, 
which is analytically solvable, has been considered by the authors of reference [10] to 
demonstrate that the RPA correlation formula works well. 

Let us define creation and annihilation operators af and a iy i = 1,2, according to : 
Xl= vb (fll + a ^ A=*^|W-ai), (2.2) 



X 2 = -L=(a 2 + a+) P 2 = i\ ^(<4 - a 2 ) ■ (2.3) 



'2SI v £ ' V 2 

The transverse frequency u, the radial frequency Q and the condensate < X 2 > will be 
determined self-consistently. Using the analogy with the linear sigma model in quantum 
field theory, the Xi and X 2 modes represent the pion and the sigma fields respectively. 
In terms of the operators dj and af, the Hamiltonian reads : 

H = Pii a t a i + Pio( a i a i + a i a i) 

i=l,2 i=l,2 

+g n (af + ai )* + g 12 (af + ai ) 2 (af + a 2 ) 2 + g 22 (af + a 2 ) 4 ' ( 2 ' 4 ) 
+h 12 {af + ai) 2 (aj + a 2 ) + h 22 (af + a 2 f + 77(0^ + a 2 ) + C 

where C is the following constant 

C = t + t+7 1 + 7?T + ? < X 2 > 2 +7^ < X 2 > 2 +^ < X 2 > 2 +p <X 2 > 4 + V <X 2 > . 
4 4 4u; 4H 2 Q a; 

(2.5) 



We have defined the following quantities : 



Pn = f + f + ^<^> 2 , (2-6) 

*> = T + & + *<*>'- <"> 

P10 = ^(pn - , P20 = ^(^22 - fi) , (2.8) 

^ n " 4^ ' " 2^ ' fe " 4^ ' (2 ' 9) 

^-^ <X2> '^ = ^ <Xa> ' (2 - 10) 



77 = 

V2tt 



= (r ] + f i<X 2 >+4g<X 2 > 3 ) . (2.11) 
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2.1 Mean Field Equations 



Normal ordering with respect to aj and af gives the result : 



H = Pii a t a i + Pio( a t a t + a i a i) 

i=l,2 i=l,2 



i=l,2 i=l,2 

+<7n { a i" 4 + a i + Qafafaiai + 4a^a^a^ai + Aafaiaiai + 12afai + Qa^af + 6a±ai j 
+#22 {°2~ 4 + a 2 + Qatata 2 a 2 + Aat at at a 2 + Aat a 2 a 2 a 2 + 12a^a 2 + Qa 2 a 2 + Qa 2 a 2 > 

, f +2 +2 , +2 2 , 2 +2 , 2 2 , o +2 + , o 2 + 

+#i2 \aj a 2 + a{ a 2 + a x a 2 + a x a 2 + 2a^ aja2 + 2a 1 a 2 a 2 
+2a^"a 1 a 2 f2 + Aafaia^a 2 + a^ 2 + a 2 + 2a ] f ~a 1 + a\ 2 + a 2 , + 20^02} 
+h\2 {a\ 2 at + a\ 2 a 2 + cl\ol 2 + a\a 2 + 2ata\at + 2a+aia 2 j 
+(/ii2 + 3/i 2 2 + ^)(a 2 + at) + Ejjfb 



(2.13) 



where 

Ehfb i + i + ^ + ^ + ( v 2 + TT + ^ + ^ <X2> J <X2> 

, 3c/ 3# c? - 

+TT^ + 777? + o - ^ + V < X 2> 
Auj 2 AO? 2ujVL 

is the mean field energy. 

Minimization of Ehfb with respect to < X 2 > gives the equation for the condensate 



77 + 3h 22 + /112 = (2.14) 

that is 

77+ < X 2 > 

Minimization of E HFB with respect to uj and Q gives the two gap equations 



H + Ag <X 2 > + — H 

\l UJ 



(2.15) 



c 2 = / i + 4^<X2> 2 +^ + ^ ! (2.16) 



n 2 = f, + 12g<X 2 > 2 +^ + ^- . (2.17) 

uj 1 Z 

We check that the coefficient of the linear term in the expression (2.12) gives the equation 
for the condensate (2.14). 

Equations (2.15)-(2.17) coincide with those written by Stevenson in [12]. When 
there is no explicit symmetry breaking (77 = 0), we have two solutions : one with < 
X 2 >= and one where the symmetry is spontaneously broken < X 2 >^ 0. For the 
symmetric solution, we have uj = Q and the gap equation is 

uj 3 -fiuj - 8g = . (2.18) 

We can compare the mean field results with exact numerical calculations [13]. For instance 
for /i — 1 and g — 1, the energies of the ground state and the two first excited states 
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are in the mean-field approximation : E = 1, 74015, E 1 = 3,90645 and E 2 = 6,07275 
whereas exact numerical calculations give : E = 1,7242,^ = 3,8304 and E 2 = 6,214. 
However, we are not aware of numerical calculations for the two-dimensional anharmonic 
oscillator in the case < X 2 >^ 0. 

2.2 RPA equations from the equation of motion method 

The symmetry generator, i. e. the angular momentum operator around the 3-axis, is given 
by: 

L 3 = X l P 2 - (X 2 + < X 2 >)P 1 , (2.19) 
or in terms of the creation and annihilation operators : 

~ ~ ' " (2.20) 

To derive the RPA equations, we will first use the equation of motion method due 
to Rowe [5]. We assume that an exact eigenstate \u > of the Hamiltonian can be created 
from the exact vacuum |0 > by an excitation operator : 

|z/>=Q+|0> and Q u \0 >= (2.21) 

Minimization of the energy E u =< u\H\u > / < v\v > with respect to a variation 5Q U of 
the operator Q v leads to the following set of equations : 

< 0| [<$<?„, [H,Q+]]\0 >= n v < 0\[5Q u ,Q+}\0 > , (2.22) 

where Q u = E u — E is the excitation energy. One has also the supplementary condition : 

<0|[ J ff,Q I/ ]|0>=0, (2.23) 

which is equivalent to generalized mean field equations [11]. 

We will restrict our choice of excitation operators to those which contain the 
same operators that appear in the symmetry generator (2.20), that is : 

Q+ = U^ ] a+ - V^cn + U [2) ata+ - V^ 2) ai a 2 + U^a^ - V^ ] a+a 2 . (2.24) 

The RPA ground state is defined by 

Q U \RPA >= , (2.25) 
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and the expectation values appearing in eqs. (2.20) and (2.23) are taken in the RPA 
ground state. RPA equations (2.20) become matrix equations which allow to determine 
the excitation energy fl u and the amplitudes U v and V v : 



A B * 
-B -A 



3— (vi 



(2.26) 



where we have defined the three-dimensional vectors : 

u v = up 

The norm matrix is given by : 

(M A M B 
\Mb Ma 

where the 3x3 matrices Ma and Mb are : 

(I \ 



U (3 v 



(2.27) 



(2.28) 



Ma = 1 + n + r 2 «i 

V «i r 1 -r 2 J 

We have introduced the notations : 

«i =< a{a{ >=< aiai > 



/0 \ 

-Ka 
\0 K 2 / 



(2.29) 



k 2 =< a 2 a 2 >=< a 2 a 2 > 



(2.30) 

ri =< afai > , r 2 =< a 2 f a 2 > (2.31) 

We have assumed that the previous quantities are real. We have also used by definition 

< a% >=< a 2 >= 0. 

The RPA basis being complete, we have the following expression for any operator 
O in terms of the RPA excitation operators and Q u : 



o= E Qt([QvM)-Qv([QtM) » 



u=l,3 



(2.32) 



where the expectation values are taken on the RPA ground state. We have therefore the 
expression of the operators af, a^a 2 and afa 2 in terms of and Q u and we deduce : 



< a~l >=< a^a^ >=< a\a 2 >= . 
We will introduce the following quantities : 

r = Ti + Ki + r 2 + K 2 + 1 , 
A = 7i + Ki - r 2 - K 2 . 



(2.33) 

(2.34) 
(2.35) 
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An important condition to remember in the EOM approach for RPA equations is 
that the norm matrix M has to be invertible in order to be able to write the normalization 
condition for the excited state \v >. 

The matrix elements of A and B are defined according to : 



An=<[a u [H,at]]> (2.36) 

A 12 =<[a 1 ,[H,af4]]> (2.37) 

A 13 =<[a 1 ,[H,a 1 4]]> (2.38) 

B u = - <[a u [H,a 1 ]]> (2.39) 

B 12 = - <[a 1 ,[H,a 1 a 2 ]]> (2.40) 

B 13 = -<[a 1 ,[H,ata 2 ]]> (2.41) 

A 2 i=<[aia 2 ,[H,at]]> (2.42) 

A 22 =< [aia 2 , [H, a+a+]] > (2.43) 

A 23 =< [aia 2 , [H, <n<4]] > (2.44) 

821 = - <[a 1 a 2 ,[H,a 1 ]]> (2.45) 

B 22 = - < [aia 2 , [if, aia 2 ]] > (2.46) 

B 23 = -< [aia 2 , [if, a+a 2 ]] > (2.47) 

^31 =< Ka 2 , [if, a+]] > (2.48) 

^ 32 =< [a+a 2 ,[H,atat\] > (2.49) 

^33 =< [afa 2 , [H,a!<4]] > (2.50) 

B31 = - < [at 02, [H, ai]] > (2.51) 

i5 32 = - < [ata 2 , [H, ai a 2 }} > (2.52) 

B 33 = - < [ata 2 , [if, a^a 2 ]] > (2.53) 



We need to calculate all these commutators. In the renormalized RPA approxi- 
mation, when calculating the expectation values in the RPA ground state, one uses the 
following approximation : 

< aiO>ja>kai >—< atiOtj >< ctkOti > + < a^ak >< OijOn > + < a^ai >< (XjOt^ > , (2.54) 



8 



where the a operators are either creation or annihilation operators of the bosons 1 and 2. 
The expectation values are calculated self-consistently whereas in standard RPA, they are 
calculated in the HFB ground state. We obtain the following expressions for the matrix 



elements of A and B : 

An = £ n = p n + 120n(r + A) + 2g 12 (T - A) (2.55) 

A 12 = fj + 3h 22 (T- A) + h 12 (A + 3T) (2.56) 

A 13 = 2h 12 A (2.57) 

A 21 = 2h 12 T (2.58) 

A 22 = (S w + £ ff )(Tl + T 2 + 1) + K lX n + K 2X a + ^12^ (2.59) 

A 23 = - S a ) - (1 + n + r 2 ) X . + % 12 rA (2.60) 

A 31 = 2h 12 A (2.61) 

A 32 = Kl (£ n + + fa - r 2 )x. + % 12 rA (2.62) 

^33 = (t 2 - Ti)(5 w - £ CT ) - K^tt - K 2 Xa + 4# 12 A 2 (2.63) 

5n = Xtt = 2pio + 2^n(r + A) + 2^ 12 (r - A) (2.64) 

B 12 = 2h 12 T (2.65) 

B\3 = V + h 12 (T + 3A) + 3h 22 (T - A) (2.66) 

B 2 i = 2/i 12 r (2.67) 

B 22 = KiXa + K 2 Xn + 4g 12 r 2 (2.68) 

B 23 = -K 2 {£ n - E c ) + (1 + n + t 2 ) Xct + 4(? 12 rA (2.69) 

#31 = 2/i 12 A (2.70) 

B 32 = -K 2 (£ n + S a ) + (n - t 2 ) X(T + 4(? 12 rA (2.71) 

S33 = ^X* + K1X0 + 4g 12 A 2 (2.72) 

The quantities E- K ,E a ,x-n and Xa are defined according to : 

£ v =<[a u [H,at]]> (2.73) 

£ CT =< [a 2 ,[#,a+]] > (2.74) 

X* = - <[ai,[H, ai ]]> (2.75) 

X. = - < [02, a 2 ]] > (2.76) 
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The equations for < X 2 > and the frequencies uj and Q (or for £ n and £ a ) are 
respectively obtained by writing the generalized mean-field equations (2.23) : 

< [H, a 2 ] >= (2.77) 

< [H, a+a+] >= (2.78) 

< [H, a+a+] >= (2.79) 

They are equivalent to the minimization of the generalized mean-field energy < H > with 
respect to < X 2 >, uo and ( see the expression of < H > eq.(2.130) below). We obtain 

fj + 3h 22 (T- A) + h 12 (T + A) = (2.80) 

PnKi +Pi (2ti + 1) + 6^11 (r + A) 2 + <7 12 (r + A)(r - A) = (2.81) 

P22K 2 + p 20 (2r 2 + 1) + g 12 (F + A)(F - A) + 6g 22 (T - A) 2 = (2.82) 

By using the definitions of Pio,P20,£w and £ a , the two last equations, which we call the 
generalized gap equations, can be written as : 

£ w (r + A) = u(2n + 1) (2.83) 

£ a (T - A) = Q(2r 2 + 1) (2.84) 

When using the generalized mean-field equations (2.80), (2.83) and (2.84), we check 
that the matrices A and B are symmetric. 



2.3 Standard RPA approximation 



In the standard RPA approximation, all expectation values are taken in the HFB ground 
state. We therefore have : T\ — r 2 — K\ — k 2 — and T = 1, A = 0. The generalized 
mean-field equations (2.80), (2.83) and (2.84) reduce to the mean-field equations given in 
subsection (2.1) , eqs. (2.15), (2.16) and (2.17). We have : £ n — cu,£ a — Q,Xn — Xa — 0. 
The matrices A and B become much simpler : 



A 



( ou 2h l2 0\ 

2h 12 u + n + Ag 12 

V o o o/ 



B 



( 

2h 12 

V o 



2h 12 

4^12 





0\ 




(2.85) 



where we have used the equation for < X 2 > : r\ + 3h 22 + h 12 

fuj 0\ 
A= oj + tt +B . 
v 0/ 



0. We note that 



(2.86) 
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The matrix M is diagonal : 





fl 









(° 





o\ 







1 


















Vo 









y0 





0/ 



(2.87) 



We see that the amplitudes and respectively in front of the operators a\a 2 and 
a\a 2 in the excitation operator (2.24) decouple. In the symmetric phase, we have 
L 3 = i(aia 2 — afa 2 ) and the operators which decouple are those which appear in the 
symmetry generator. 

The standard RPA equations therefore reduce to a 4 x 4 system (instead of 6 by 6). 
Its dimension will be then divided by two and the RPA frequencies fl u satisfy : 

det [(A -B)(A + B)- VtJ] = , (2.88) 

where now A and B are 2 by 2 matrices and I is the 2 by 2 unity matrix. This gives : 



(c 2 - nl 



CO 



+ 8g 12 (uj + n))- i6/i 2 2 ^ + n) = o 



(2.89) 



This equation is valid for the two solutions : < X 2 >= and < X 2 >^ 0. 

Let us check the existence of a zero RPA frequency in the case of a ground state 
with broken symmetry (the analog of the Goldstone mode in the linear sigma model). (In 
the following, we have no explicit symmetry breaking : 77 = 0). By introducing the two 
following quantities (which correspond to loop integrals in quantum field theories) : 

h = %u~ ' Ia = 2n ' 

the two gap equations in the case < X 2 >^ can be written as : 



(2.90) 



(2.91) 

n 2 = 8g< X 2 > 2 . (2.92) 

We then introduce the following quantity (which is analog to the self-energy in quantum 
field theories) : 

ml) = " + ■ (2.93) 



2luH ttl-(cu + tt) 2 ' 
Equation (2.89) for the RPA frequencies can be rewritten as : 

u + Q 1 



(Ql-u 2 ) (l-4g L - 



- 32g 2 < X 2 > 2 



or 



Ug 2 < X 2 > 2 



8gmi) 



= 0, 

(2.94) 
(2.95) 
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This equation is identical to eq. (28) of [6]. We also notice that : 

s(n " = °> = -2^Tnr (2 ' 96) 

h ~ l " = %Mu + a) ■ (2 - 97) 

Therefore : 

E(n, = 0) = . (2.98) 

We then use the gap equation (2.92) to write the equation for the RPA frequencies (2.95) 
as : 

and then we use the first gap equation (2.91) to obtain : 

^ = x _^ mi) ^ 2 ~ (ml) - ml = o)) • (2.100) 

It is clear now that there is a zero frequency in the RPA spectrum. 



2.4 Renormalized RPA approximation 



We see from eqs. (2.29) and (2.30) that if k x or k 2 are not vanishing, the norm matrix A/" 
is not diagonal. In this case, it won't always be possible to write normalization conditions 
for the excited states \v >. We therefore impose the supplementary conditions : < 
afaf >=< a±ai >= and < a^a^ >=< a 2 a 2 >= 0, that is k± — and k 2 = 0, i.e. 
we neglect pair correlations. With these conditions, the matrix Af stays diagonal in the 
renormalized approximation : 




with Ma 



o 





T\ 





T2 






Ti - r 2 J 



The matrices A and B simplify : 
A= 







(w + n)(l + ti + r 2 ) 









+ B 



B = 



( ° 

2h 12 r 

\2h 12 A 



2h l2 T 



%i 2 r 2 

Ag 12 TA 



(n-w)(ri -T 2 )j 

2h 12 A \ 
4g 12 TA . 
Ag 12 A 2 J 



(2.101) 



(2.102) 



(2.103) 
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We notice that, contrary to what happen in standard RPA, the operators a\a 2 and afa 2 
don't decouple : the RPA matrix remains 6 by 6. Its expression is very similar to standard 
RPA at finite temperature, T\ and r 2 being the occupation numbers [9]. 

We defined new three-dimensional vectors U u and V v and a new RPA matrix 7Z by 



'U v 



M 1 ' 2 ^ , ^ = ^- 1/2 ( K A B B A )- N " 1/2 - ( 2 - 104 ) 



We have : 



K 



'A B N 
-B -A] 

with : A = M A 1/2 AM A 1/2 and B = Ma~ 1/2 BMa~ 1/2 ■ 
then write : 



A B 
[-B -A j 

and the RPA frequencies are determined by : 



Vj \V V 



(2.105) 

The renormalized RPA equations 
f U v \ 



det 



(A -b)(a + b)- n 2 j] = , 



I being the 3 by 3 unity matrix. The matrices A and B are given by : 

(u 0\ 



A 



B = 



w + fl 

\0 Sl-wJ 

( 2h 12 VT 2h 12 VA \ 

2h 12 VT 4g 12 T 4g 12 VTA 

\2h 12 y/A 4g l2 VTA 4g 12 A ) 
where r = 1 + T\ + r 2 and A = T\ — r 2 . 

The norm of the excited states 



< u \ u >=< 0\Q u Q+\0 >=< 0\[Q v ,Q+]\0 > , 
can be chosen to be equal to one. This corresponds to : 



or 



i.e. 



M A 














(0 





( t (i \ ni.\ 

m 2 + m 2) r + m 2 - m 2 - m 2 - ^ 



(2.106) 



(2.107) 



(2.108) 



(2.109) 



(2.110) 

(2.111) 

(2.112) 
(2.113) 
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By replacing the expressions for g± 2 and h± 2 , equation (2.107) can be written in the 



form : 



4g 



uj + Q 



uj — Q 



-32g 2 < X 2 > 2 



(u + ay 

uj — Q 



con ni-(to-n)\ 

A UJ + tt 

+ 



r 



We introduce the quantity : 

2N uj + Q 1 + ti + r 2 



UJ — Q T\ — T 2 



(2.114) 



(2.115) 



2ujQ Q 2 u -(uj + Q) 2 2ujQ Q 2 - {uj - Q) 2 ' 

which is formally similar to the self-energy operator in quantum field theories at finite 
temperature. Equation (2.114) then becomes : 

Q 2 u = uj 2 + 64g 2 < X 2 > 2 



1 (2 - U6) 

This equation for the RPA frequencies VL V has the same form as in standard RPA (eq. 
(2.95)) but, in renormalized RPA, S r contains the densities T\ and r 2 which have to be 
determined self-consistently. 

Proceeding in the same way as in standard RPA, let us check the existence of a 
zero frequency in renormalized RPA in the case < X 2 >^ 0. We first introduce I w = 
(2ti + l)/2u and I a = (2r 2 + l)/2fi. The equation for the condensate and the two 
generalized gap equations have then the same form as in standard RPA : 



and we have again 



We have : 



or 



E r (nl = o) 



H + Ag < X 2 > z +Agh + I2gl a = 

uj 2 = ii + Ag < X 2 > 2 +12gl n + 4gl a 
n 2 = /2+12g< X 2 > 2 +AgI 7T + \2I a 

uj 2 = 8g(I n - I a ) 
VI 2 = 8g< X 2 > 2 . 

l l + n + T 2 

~2ujVL (oj + Q) 

i tt In 



l n - t 2 

2con (uj - ft) 



MK = o) 



UJ 



ft 2 



(2.117) 

(2.118) 
(2.119) 

(2.120) 
(2.121) 

(2.122) 

(2.123) 



Using the two equations (2.120) and (2.121), we obtain the same expression as in standard 
RPA with E replaced by S r : 

8g 



ft 



(n 2 -uj 2 ) (E r (n2) - E r (rs = o)) 
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(2.124) 



We have therefore proven the existence of a zero frequency in the RPA spectrum in 
the case where we have a broken symmetry in standard RPA as well as in renormalized 
RPA. This result is very encouraging to apply renormalized RPA in the linear sigma 
model. 

The operator associated to the zero mode is hermitian and can not be normalized 
according to (2.113). It is equal to the symmetry generator L 3 . The amplitudes U and 
V for the zero mode are : 

(2.125) 
(2.126) 

C^h|- + ^]VA, (2.127) 
and Uq = Vo- 

In order to obtain a closing of the renormalized RPA eigenvalues problem, one needs 
the expressions of the expectation values T\ and r 2 in terms of the RPA amplitudes. One 
uses the inversion formula (2.32) and the algebra of the sp(4) group : 




aJ ai 
1 2 



^ ([aia 2 , afa%] + [afa 2 , a^] - l) , (2.128) 

a t a 2 — 2 (t a i a 2) at a t] ~ [ a i a 2, Q>iQ>t] ~ l) • (2.129) 

This provides a set of supplementary equations and allows a closing of the renormalized 
RPA eigenvalues problem. 

To be complete, we finally give the expression of the energy in the renormalized 
RPA approximation. From the expression of H with normal ordered product (2.12), we 
obtain : 

< H >= Pn(Ti + «i) + p 22 (r 2 + K 2 ) - UJK\ ~ 

+ 12g U (T 1 + K 1 )(T 1 + K 1 + l) + 12g 22 (T 2 + K 2 )(T 2 + K 2 + l) . (2.130) 

+2g 12 (2(r 1 + «i)(r 2 + k 2 ) +t 1 + k 1 + t 2 + k 2 ) + E HB 

We have checked that minimization of < H > with respect to < X 2 >, uj and Q gives the 
equation (2.80) for < X 2 > and the two generalized gap equations (2.81) and (2.82). 

2.5 Renormalized RPA for the symmetric solution 

For the symmetric solution < X 2 >= 0, we have u = f2, T\ = r 2 . We still use K\ — k 2 — 0. 
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The matrix Af is diagonal with a vanishing matrix Mb and : 

(I 0\ 
1 + 2n 



The renormalised RPA matrix writes : 

/uj 0\ 
.A = 2w 
v y 







(2.131) 



/o 



+ B with £ 







0\ 



o Sa + ^i) 







(2.132) 



For the symmetric solution we therefore obtain a decoupling of the amplitudes in front 
of the operators aia^ and afa2 in renormalized RPA, contrary to what happens for the 
solution with broken symmetry. 



— (i + 2n) 



The RPA frequencies are the solutions of : 

(a, - uj 2 ) 

and the generalized gap equation writes : 

= /i + — (2n + l) 



= 



(2.133) 



(2.134) 



We have therefore the RPA frequencies : Vl 2 = uo 2 and fi 2 , = 5cj 2 — //. 



3 RPA from the time-dependent formalism 

In this second part of our paper, we will derive the RPA frequencies from the linearization 
of the time-dependent Hartree-Bogoliubov (TDHB) equations. This approach has been 
introduced in many-body non-relativistic theories [14]. It has also been used in A$ 4 field 
theory in references [15], where small oscillations in the broken phase lead to one and 
two meson modes of the theory. In this formalism, the natural variables are < X i > and 
< Pi > (or $(f) and Tl(x) in A$ 4 field theory [15]). This formalism is well adapted to 
dynamical problems. Some people working on RPA approximations use the formalism 
with the creation and annihilation operators, other use the time-dependent variational 
approach. It is interesting to make a close comparison between the two approaches. 

For a two-dimensional system, a Gaussian state at finite temperature can be de- 
scribed by a vector a a and a matrix E ab , a — 1,2: 

a _( x a \ _ ab _ ( 2G ab -iT ab \ 

\ -ip a J ' " V ~iT ba -2S ab ) ^ ' ' 
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with 

-a =<X a > , p a =< P a > (3.2) 

Q ab =< xaj^b > (33) 
S ab =< papb > (3^ 
T ab =< jfapb + pb^a > (3 ^ 

with X a = X a - < X a > and P a = P a - < P a >. The matrices G and S are symmetric. 

At zero temperature, we have only two independent matrices among G,T and S. 
The state can be described by a Gaussian wave function parameterized by x, p, G and E : 

ij(X\X 2 ,t) = ^exp(^- <X -x\— + i^X -x>^exp(i<p\X -x>) , (3.6) 

where the matrix E is related to the preceeding matrices by : 

T ab = 2(GE + m) ab (3.7) 

S ab = ^(G- 1 )^ + 4(EGE) afe (3.8) 

If we work with the operators ai,af, as it is more usual in many-body problems, 
one introduces the matrix p defined by : 

l + 2p ) lJ = 3 , I . 3 ) 3.9 

\ 2 < ajaj > - < ajdj + CLjaJ > J 

The link with the representation with operators Oj, af and the representation with Xi, Pi 
is given by the matrix relation : 

S' = i(l + r)(l + 2p)(l + T) (3.10) 

where 



I being the 2 by 2 matrix and 



< XiXj > ^uJ- -iy/cJl < XiPj + PjXi > J 
-iJul < X.Pi + PiX, > 4= -24= < PiPj > 4= 

with uji = uj and u 2 = £1 (see eqs. (2. 2), (2. 3)). 



(3.12) 
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3.1 The time-dependent Hartree-Bogoliubov equations 



In this part of the paper, we will use the following notation for the Hamiltonian of two- 
dimensional quantum anharmonic oscillator : 

H = &&^ + x !-^ 2 - < 313 > 

The time-dependent variational Hartree-Bogoliubov equations are obtained at zero 
temperature from the minimization of : 

S = J dt < ip\id t - H\ip > , (3.14) 

where \ip > is the state corresponding to the variational wave function (3.6). At finite 
temperature, the TDHB equations are obtained by minimizing [16] 

Z(V(t)) = triVfa)) - jT 1 * tr + i[H, V(t)]j , (3.15) 

where T>(t) is a variational density matrix, chosen to be a Gaussian and therefore char- 
acterized by a and S. 

The TDHB equations can be written in the following compact form [17] : 

ia = tw , (3.16) 

<3 = -[(S + r)W(S-r)-(S-r)W(S + r)] , (3.17) 

or 

iE = 2[EHr -tHE] . (3.18) 
The vector w and the matrix Ti are defined by : 

6 < H >= w?5a? - hr (HpZ%) . (3.19) 

For the Hamiltonian (3.13), we have : 

1 



<H>= ^-(pj+pl + trS) 



+gG u (6xl + 2x 2 2 - 2a 2 + 3G n + G 22 ) 
+gG 22 (2x\ + Qxj - 2a 2 + G u + 3G 22 ) 
+^G 12 (8^ 1 S 2 + 4G 12 ) + g(xj + x\- a 2 ) 2 



(3.20) 



The vector w and the matrix TC are given by : 

i 



w 2 = —pi , w 2 = —p 2 (3.21) 
m m 
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w 



w 



1 = Agxt (x\ + x 2 2 - a 2 ) + 12g Xl G n + Ag Xl G 22 + 8gx 2 G 12 (3.22) 

2 = Agx 2 (x\ + x\- a 2 ) + \2gx 2 G 22 + Agx 2 G 22 + Sg^G 12 (3.23) 

n\l = ~6gx 2 i - - 6^G n - 2gG 22 + 2ga 2 (3.24) 
Ull = -4gxix 2 - AgG 12 = n 2 l (3-25) 
nH = ~6gxl - 2gx\ - QgG 22 - 2gG u + 2ga 2 (3.26) 

s ^ = ^ (3 - 27) 

«;5-2i^^ = , W g = 2*^-^ = (3.28) 

3.2 Static solution of the TDHB equations 

The static solution of the TDHB equations (3.16) and (3.18) is given by : 

w = (3.29) 

En(a,^)T-rn(a,E)E = (3.30) 

Using the rotational invariance we can choose x 1 = and from w = we deduce also 
G 12 = (to simplify the notations we don't use the bar on the matrix elements of S for 
the HFB ground state). From w = 0, we deduce the existence of two solutions. One 
solution is symmetric with x 2 = 0. The other solution shows a broken symmetry with 
x 2 7^ and given by : 

xl = a 2 - G 11 - 3G 22 (3.31) 

/,From eq. (3.30), we deduce : 

T ab = (3.32) 
S 12 = (3.33) 
S 11 = AmgG 11 {x\ - a 2 + 3G 11 + G 22 ) (3.34) 
S 22 = AmgG 22 (3x 2 2 - a 2 + G n + G 22 ) (3.35) 
For the symmetric solution, we have : 

For the solution with broken symmetry, 

S 11 = 8gm(G u - G 22 )G 11 (3.37) 
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S 22 = 8gmG 22 x 2 2 (3.38) 

At zero temperature, we have for the static solution : T ab = and S 11 = 1 /AG 11 , S 22 = 
1/AG 22 . Equations (3.34) and (3.35) give the two gap equations : 

^Gn = 2g(x 2 2 - a 2 + 3G n + G 22 ) (3.39) 

^G 2 2 = 2g(Zx 2 2 -a 2 + G u + G 22 ) (3.40) 

By using the identification G 11 = l/2u and G 22 = 1/2S1, the previous equations are 
identical to the mean- field gap equations obtained in the first part of the paper (2.16) 
and (2.17). 

At zero temperature, the symmetric phase is characterized by : 

x 2 = (3.41) 

G u = G 22 , ^G u 2 = 2g(AG 11 - a 2 ) . (3.42) 
The solution with broken symmetry at zero temperature is characterized by : 

x\ = a 2 - G n - 3G 22 (3.43) 
! G u 2 = g(G 11 - G 22 ) (3.44) 



32m 



' G' 2 2 = gx 2 2 (3.45) 



32m 



3.3 Small oscillations around the static solution 

The linearization of the TDHB equations (3.16) and (3.18) writes 

i5a = 5w (3.46) 

i5E = 2 [8E H t - t H 5E + E 5H r - r 5H s] (3.47) 

where H is the matrix Ti evaluated for the HFB static solution d,E. 

Let us write more explicitly the linearization of the TDHB equations around the 
static HFB solution with broken symmetry characterized by x\ = 0, x 2 ^ 0, G 11 , G 22 , S 11 , S 22 ,pi 
p 2 = 0,G 12 = S 12 = T ab = 0. We obtain a differential system with 14 variables : 
5x u 5 Pl , 5x 2 , 5p 2 , 5G 11 , 5G 12 , 5G 22 , 5T 11 , 5T 12 , 5T 21 , 5T 22 , 5S 11 , 5S 12 , 5S 22 : 

8x 1 = —5p! (3.48) 
m 
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8p x = -8g(G n - G 22 )5 Xl - 8gx 2 5G 12 (3.49) 

5x 2 = — 8p 2 (3.50) 
m 

Sp 2 = Sgx 2 2 - Ylgx 2 bG 22 (3.51) 

SG U = —5T U (3.52) 
m 

SG 12 = —5T 12 + —5T 21 (3.53) 

2m 2m v ' 

SG 22 = 1st 22 (3.54) 
m 

5T n = -lQg(G u - G 22 )5G n + -5S 11 (3.55) 

m 

5f 12 = -l6gG n x 2 5x 1 - 16g(x 2 2 + G n )5G 12 + -5S 12 (3.56) 

m 

5f 21 = -16gG 22 x 2 6x 1 - lQgG u 5G 12 + -5S 12 (3.57) 

m 

5f 22 = -16gx 2 2 5G 22 + -5S 22 (3.58) 

m 

5S n = -8g(G n - G 22 )5T n (3.59) 

SS 12 = -4g(G u - G 22 )5T 12 - 4gx 2 2 5T 21 (3.60) 

SS 22 = -8gx 2 2 5T 22 (3.61) 



To write these equations, we have used the mean-field equations for the solution with 
broken symmetry. 

^From this differential system of first order, we obtain a 14 by 14 RPA matrix. 
Coming back to the operator representation, with the help of relation (3.10) and its 
inverse, this corresponds to the 14 operators we have to include in the most general form of 
the excitation operator of the first part of this paper if we keep only bilinear operators : 
ai, af, a 2 , a 2 , a±ai, afaf, a 2 a 2 , a 2 a 2 , a±a 2 , afa 2 , afa 2 , a\a 2 and afa±, a 2 a 2 . The last two 
operators appear only at finite temperature (For the derivation of RPA equations at finite 
temperature, see reference [18]). 

The operator L 3 allows to separate these operators in two sectors : the six op- 
erators corresponding to the "pion" sector and the 8 op- 
erators a 2 , a 2 , aiai, afaf, a 2 a 2 , a 2 a 2 , a\a 2) afa 2 , afa 2 , aia 2 , afai, a 2 a 2 corresponding to 
the "sigma" sector. By using again the correspondence between the Xi,Pi representation 
and the af,di representation, we deduce that the variables in our first-order differential 
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system which correspond to the "pion" are : 5x 1 ,5pi,5G 12 ,5T 12 ,5T 21 ,5S 12 , the other 
corresponding to the "sigma" sector. We check indeed that the RPA matrix, which we 
call 1Z, can be written in two blocks, one 6 by 6 corresponding to the "pion" and one 8 
by 8 corresponding to the "sigma" , the two sectors being disconnected. We can write the 
first-order differential system in the form : 

5X = 7Z5X , (3.62) 

where 5X = (6x u Sp u 5G 12 , <5T 12 , 5T 21 , 5S 12 , 5x 2 , 5p 2 , 5G 11 , 5G 22 , 5T U , 5T 22 , 5S 11 , 5S 22 ) and 

V o mJ 

The 6 by 6 RPA matrix M. n for the "pion" sector is equal to : 



(3.63) 



M 7 



M a = 



( 





J_ 

m 











o\ 


-8g(G 


11 _ 


G 22 ) 


-8gx 2 






















i 

2m 


i 

2m 





-16gG u 


x 2 


-16g(3* + G 11 ) 








2^ 

m 


-16gG 22 


x 2 


-16^G n 








2_ 

m 


\ 











-Ag(G u - G 22 ) 


-4gx 2 


0/ 


3 by 8 RPA matrix M a 


is equal to : 








( ° 


l_ 

m 

















-8gx 2 








-12^2 























j_ 

m 























J_ 

m 











-16^(G n 


- G 22 ) 








_2_ 

m 











-16^2 























-8g(G n - G 22 ) 








V o 














-8gx 2 2 






(3.64) 



0\ 









_2_ 

m 





0/ 

(3.65) 

We have detA^ = 0. However, before to conclude about the existence of a zero mode 
associated to the spontaneously breakdown of the rotational symmetry, we have to elim- 
inate the spurious modes corresponding to invariants of the TDHB evolution. At zero 
temperature, this will allow to reduce the M w matrix to a 4 by 4 matrix. This is also the 
dimension we have found in the first section in standard RPA. 



At zero temperature, we have the following condition to be satisfied to have a pure 



state 



or 



p(p + /)=0 



-S'rS' = T 



(3.66) 
(3.67) 
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where the matrices p and S' are given by eqs. (3.9) and (3.12) and I is the 2 by 2 unity 
matrix. We linearize eq. (3.67) around the HFB static solution. This gives the following 
conditions : 

G ik 5T jk - 5T ik G kj = (3.68) 

k 

ST ki S kj - S ik 5T kj = (3.69) 

k 

G tk 5S kj + 5G ik S k3 = (3.70) 

k 

Equations (3.68) and (3.69) give respectively : 

G n ST 21 - 5T l2 G 22 = (3.71) 

S 22 5T 21 - 5T l2 S n = (3.72) 

At the minimum at zero temperature, we have : S* 11 = 1/AG 11 and S 22 = 1 /AG 22 . The 
two previous conditions are therefore equivalent. ^From eq. (3.70), we obtain : 

5G 12 + AG n G 22 5S 12 = (3.73) 

5G n +A(G n ) 2 5S n = (3.74) 

5G 22 + 4(G 22 ) 2 5S 22 = (3.75) 

The first condition is for the "pion" sector and the last two conditions are for the 
"sigma" sector. 

Coming back to the operator representation, we have : 

. _ _ 2 ~ ~ 

2Vutt < XxX 2 > < P X P 2 >= 2(< ai a% > + < afa 2 >) (3.76) 

< XxPs + P 2 X 1 > -2d- < X 2 Pi + PiX 2 >= 4(< aia+ > - < a+a 2 >) (3.77) 
V s i V u 

By using G 12 =< X^ 2 >,S 12 =< P X P 2 >,T 12 = 2 < X X P 2 >,T 21 = 2 < X 2 P 1 >, and 

G 11 = l/a;,G 22 = I/O, for the static HFB solution, we see that conditions (3.71) and 

(3.73) correspond to the decoupling of the operators aia 2 , afa 2 we have found in the first 

section in standard RPA. 

For the pion sector, we will therefore consider the following new variables : 

SY = ( 6x u 5 Pu 5C = -AG n G 22 5G 12 + 5S 12 , 5D = G U 5T 12 + G 22 5T 21 , 

5E = 5G 12 + G 22 5T 12 - G U 5T 21 + 4G n G 22 5S 12 , (3.78) 
SF = 5G 12 - G 22 5T 12 + G U 5T 21 + 4G U G 22 5S 12 ) 

We have : 

SY = M' w 5Y (3.79) 
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/,From the condition to remain in the variational space of pure Gaussian states we have 
: SE = 0, SF = 0. We therefore consider in M.' n the 4 by 4 matrix corresponding to the 
four coordinates Sxi, 5p±, 5C, SD (we notice that there are non- vanishing matrix elements 
in the fifth and sixth columns of M'^) : 



m: 



( o 

-8g(G u - G 22 ) 


V-16((G 11 ) 2 + (G 22 ) 2 )( ?a : 2 







n 32G 11 G 2I gx 2 
U l+16(G 1:l G 22 ) 2 



N 2 



\ 


/ 



(3.80) 



where 
N 1 = 

No = 



m((G u ) 2 + (G 22 ) 2 ) 
2 

"m(l + (G n G 22 ) 2 ) 



G ^G 22 {G 22 - 2gm) + (G ii ) 2 (G 22 + 2gm) + 2G 22 gx 2 2 m 



(3.81) 



G + G + 32(G LL YG 22 gm + 32G G gm(G + x 2 2 



(3.82) 

We then calculate det Ai" and, by using the two gap equations (3.44) and (3.45) for the 
solution with broken symmetry, we obtain : 



det AC = 



(3.83) 



At zero temperature, we have therefore checked the existence of a zero frequency for the 
solution with broken symmetry in the standard RPA. 



At finite temperature, the condition (3.66) becomes : 



p(p + I) = -(C-l)I 



(3.84) 



where C is called the Heisenberg invariant and is a quantity conserved by the TDHB 
evolution : C — 0. In our 0(2) model, C has two indices : C ab , a, b = 1, 2. The Heisenberg 
invariant is equal to : 



C ab = 4 < X a X c >< P c P b > - < X a P c + P c X a >< X c P b + P b X c > (3.85) 

c 

For the static solution at finite temperature, we have T ab = and : 

AG u S n = C n , 4G 22 1 S 22 = C 22 . (3.86) 

C 11 and C 22 are related to the occupation numbers for the bosons 1 and 2 according to : 
2n a + 1 = \JC aa . Similarly to the zero temperature case, we linearize the condition (3.84) 
around the static solution characterized by G 11 , G 22 , S 11 and S* 22 . It is then convenient to 
introduce new variables in the pion sector and we check the existence of the zero mode. 
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4 Conclusion 



In this paper, we have applied the random phase approximation and its extension called 
renormalized RPA to the quantum anharmonic oscillator with an 0(2) symmetry. The 
expression for the RPA matrix in renormalized RPA is formally very similar to the RPA 
matrix appearing in standard RPA at finite temperature. We focused on the existence of 
a zero mode among the RPA frequencies in the case where the ground state has a broken 
symmetry. This result is encouraging to apply renormalized RPA in the linear sigma 
model. We have compared also the approach with the creation and annihilation operators 
with the time-dependent approach and we identify the variables corresponding to the 
"pion" sector and those corresponding to the "sigma" sector. The numerical resolution 
of the self-consistent renormalized RPA equations will be the subject of a next paper, 
where Hartree-Bogoliubov mean-field results, standard RPA results and renormalized 
RPA results will be compared to exact numerical results for the vacuum energy and the 
energy of the first excited states in the cases of the symmetric solution and the solution 
with broken symmetry. 
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